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Abstract 

Although BCS pairs of fermions are known not to obey Bose-Einstein (BE) commutation relations 
nor BE statistics, we show how Cooper pairs (CPs), whether the simple original ones or the CPs recently 
generalized in a many-body Bethe-Salpeter approach, being clearly distinct from BCS pairs at least 
obey BE statistics. Hence, contrary to widespread popular belief, CPs can undergo BE condensation to 
account for superconductivity if charged, as well as for neutral-atom fermion superfluidity where CPs, 
but uncharged, are also expected to form. 
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1 Introduction 



A recent electronic analog of the Hanbury Brown- Twiss photon-effect experiment Samuelsson and Biittiker 
[T] strongly suggest electron pairs in a normal/superconductor junction to be bosons, although further work 
seems needed to provide compelling empirical proof. On the other hand, assuming Cooper pairs (CPs) 
to be bosons, it was recently proved [H [3] within a generahzed Bose-Einstein condensation (BEC) model 
that includes the BCS-Bose crossover theory as a special case, that a BCS condensate is precisely 

a Bose-Einstein (BE) condensate consisting of equal numbers of particle- and hole-CPs, at least for the 
Cooper/BCS [inillT] model interfermion interaction in the limit of weak couphng when the crossover picture 
reduces to BCS theory. This is significant since BCS argued (Ref. [T7j, footnote 18) that "owr transition 
is not analogous to a Bose-Einstein condensation (BEC)." Somewhat later, Bardeen wrote [18j that "...the 
picture by Schafroth (1955). ..of electron pairs... which at low temperature undergo a BEC, is not valid.^^ In 
contrast, both Feynman [12] and Josephson [2D] referred without proof to CPs as bosons. 



On the other hand, a long-standing common objection to attempts [51 [51 HH [H] to unify BCS and BEG 
theories for a description of superconductivity has been that CPs, strictly speaking, cannot be considered 
bosons. In this Letter we establish a clear distinction between BCS pairs and CPs, showing that while 
the former are not bosons the latter very definitely are, at least insofar as they obey BE statistics. The 
distinction, though elementary, is hardly if ever mentioned in the superconductivity literature, though it 
appears not to be as unacceptable among the neutral-fermion superfluidity community where it is also 
applied, albeit with a different interfermion interaction, to the BCS and Cooper pairs as in liquid ^He [23] 
and in ultracold trapped alkali Fermi gases such as '*°K [21] and ^Li [55]. Recent studies [2B]-[3D| have also 
appeared addressing superfluidity in degenerate Fermi gases. Indeed, a BE-like condensate of dimers formed 
from neutral '^''K [21] or ^Li [35]- [31] fermionic atoms via magnetically-tunable "Feshbach resonances" has 
only recently finally been observed. 



2 Cooper pairs 

The original or ordinary CPs [16] emerge by solving a Schrodinger-like equation in momentum space for two 
particles above the ideal-Fermi-gas (IFG) sea of the remaining TV — 2 system fermions. If two-hole pairs 
are also included, the extended Cooper eigenvalue equation yields [35] [36] purely imaginary values. This 
precludes defining CPs relative to an IFG sea. However, a more both general and self-consistent treatment 
starts with an unperturbed Hamiltonian representing not the Fermi sea but BCS-correlated ground state 
in the Bethe-Salpeter (BS) integral equations with hole propagation. When solved [37|-[39] in the ladder 
approximation for bound two-particle and two-hole states, non-purely-imaginary eigenvalues are restored for 
what has been called the "moving CP" solution of the BS equations. The BS results in 3D [37], 2D [38^, and 
ID [39] are consistent with each other, as expected. 

Unlike the Fermi-atom case, the interfermionic interactions in (especially cuprate) superconductors are a 
deep, problematic and highly controversial issue [40j-|43j. So, we confine ourselves to an electron gas inter- 
acting pairwise via the ingeniously simple Cooper/BCS model interfermion interaction jl61 117] 

y^^,^l~V iffcj. or y]4^< |k±iK|, |k'±iK| < (1) 
[ otherwise, 

where Vkk' = L^"^ J dr J dr'e~^^''^V{r,r')e'^^ is the double Fourier transform of V(r,r'), the (possibly 
nonlocal) interaction in cZ-dimensional coordinate space, with r the relative coordinate of the two electrons. 
Here V > 0, and Hlud = ?l^fc|)/2m is the maximum energy of a phonon associated with the vibrating ionic 
lattice underlying the electron gas, while Ep = h^kp/2m is the Fermi energy and m is the effective electron 
mass. In Ref. [16] there are no hole pairs (viz., below Ep) so the lower limit in ((T]) was kp, while in Ref. 
[17] it is \/kp — k^- In the simpler Cooper case it means that two electrons with momentum wavevectors 
ki and k2 interact with a constant net attraction —V when the tip of their relative-momentum wavevector 
k = §(kj^ — k2) points anywhere inside the shaded overlap volume in fc-space of the two spherical shells in 
Fig. 1 whose center-to-center separation is K = ki -|-k2, the total (or center-of-mass) momentum wavevector 
for the pair. Otherwise, there is no interaction and hence no pairing. For neutral-fermion systems there is 
no cutoff hojD and V is not constant in k and k', except when the interaction range is much shorter than 
the average fermion spacing as in a trapped Fermi gas where it has become customary to employ contact 
(or 5) interactions. Such are a special case of the more general finite-ranged, separable interaction, as is ([T]), 
introduced in Ref. [8], and for which the conclusions below will continue to apply. 

Vectors k ending at all points of a simple-cubic lattice in fc-space (in the 3D case) with lattice spacing 2tt/L 
with L the system size, ensure that the interaction ([T]) is nonzero provided such points are in the shaded 
overlap volume of Figure 1. The endpoints of two vectors k and k' are illustrated in Figure 2. In the 
thermodynamic limit there are infinitely many acceptable k values allowed for each fixed value of K for the 
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Figure 1: Cross-section of overlap volume in fc-space (shading) where the tip of the relative pair wavevector 
k = ^(ki — k2) must point for the attractive Cooper/BCS model interaction ([1]) to be nonzero for a pair of 
total (or center-of-mass) momentum KK = ?i(kj^ + k2). The tail of k is at the midpoint of K. 



K >0 CP eigenvalue equation [l6][44] with interaction H]), namely 



(2) 



The prime on the summation sign signifies the restrictions in ([T]) using kp as lower limit. For CPs consisting 
of equal and opposite momentum electrons K = 0, solving ^ analytically gives the familiar result of 
infinite-lifetime, negatively-bound (with respect to 2Ep) CPs 



So = ~2hLUD/[exp{'2/N{EF)V} - 1] 



-2/A 



(3) 



where A = N{Ef)V > with N{Ef) the density of electronic states at the Fermi energy. The first expression 
for £o is exact for all values A in 2D [where iV(e) is constant] and is otherwise a good approximation if 
hujD ^ Ef as occurs in metals. If hole CPs are included in the BS scheme along with electron CPs as in ([2]) 
and ([n]), instead of ([3]) one obtains the evidently unphysical purely- imaginary result 



So = ±i2hujDl\/e'^/>- - 1 



(4) 



(see detailed treatment in Ref. [22] )■ Replacing the IFG sea by the BCS ground-state sea, instead of ([2]) the 
eigenvalue equation for Sk in 2D [38] (and a very similar one in 3D [Z7\ and in ID [39]) is rather 



— XhvF / dk d<y5WK/2+kWK/2-k X {uK/2-kUK/2+k " WK/2+kWK/2-k 

Jkp-ko JO 



E, 



K/2+k 



E- 



K/2-k 



(5) 



where vf 
Ek 



^^K + (-E'K/2+k + -EK/2-k)^ 

\j2EFlm, if is the angle between K and k, k^ is defined by k^/kp = huj£,/2EF <C 1, 
\/^fc^ + A^ with ^fe = h^k'^/2m — Ef and A is the electronic gap, while vl = ^(1 — S,k/Ek) and 



3 



27I/L 




Figure 2: Two distinct relative- momentum wavevector values k = ^(ki — k2) and k' = ^(k'j^ — k^ (dashed 
lines) whose tips are within the overlap volume of Figure 1. They contribute to the summation in ([2|) or to 
the integral in ([5]) to determine a CP energy £k- This illustrates why the (5k, k' term for BCS pairs in ([15]) 
is always zero for CPs, with the Pauli Principle being strictly satisfied. 

u\ = \ — vl. are the usual Bogoliubov functions [IS]- Note that the lower limit in ([T]) is being taken as 
\/kp — kj) as in BCS theory [17] which does not exclude hole pairings. Besides endowing CPs with the 
physically expected finite lifetimes, the more general BS treatment yields the nontrivial eigenvalue solution 
for K = given by £o — 2A, with A = Hujd/ smh.{l/\) — > 2hujDe~^^'^ the single-electron BCS energy 

A — >0 

gap, as opposed to the ordinary CP problem giving ([3]). [Note the factor of two difference in the exponential 
compared with ([3]) due to excluding [16] or not [17] hole pairs.] The nontrivial (new) BS solution governed by 
([5]) has been called [37]- [39] a "moving CP" composite-boson excitation mode to distinguish it from another 
(a trivial or known) sound wave solution sometimes called the Anderson-Bogoliubov-Higgs ([35] p. 44), [46] . 
[47] excitation mode which in contrast with the moving CP mode vanishes as K ^ for any finite coupling, 
i.e., is gapless, and is governed by another eigenvalue equation [37] [35] different from ([S]). In the pure boson 
gas, on the other hand, both "particle" and "sound" solutions are indistinguishable 05] In the fermion 
case of interest here it should be feasible to search for, identify and distinguish both particle and sound 
modes experimentally, e.g., see Ref. [50] . 



3 BCS pairs with nonzero total momentum 

Whether the pairwise interfermion interaction is between charge carriers or between neutral atoms, a CP 
state of energy £k, defined via eigenvalue equation ([2]) or ([5]), will clearly be characterized only by a definite 
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K but not definite k. Although elementary, this is our main point. It contrasts with that of a "BCS pair" 
defined as a dimer with fixed K and k (or equivalently fixed ki and k2), even though only the case K = is 
considered in Ref. [17j . for which their annihilation 6k and creation operators are not quite bosonic since 
they obey the relations, Ref. [H] Eqs. (2.11) to (2.13), 



= (1 - n_ki - nk^)S] 



kk' 



(6) 



(7) 



where n^^^ — a^[^^a±ks are fermion number operators, with creation a]^^^. and annihilation Okis operators 
referring to the fcrmions, and 

{6k,6k'} = 26k5k'(l-<5kk') (8) 

which is not quite fermionic, unless k = k' when (jlSp is not bosonic. The precise Bose commutation relations 
are, of course, 



6k, 



- [bk, 6k'] = 



(9) 
(10) 



with ^ differing sharply from ([5]). The fermion creation aj^^^ and annihilation akjs operators satisfy the 
usual Fermi anti-commutation relations 



{^kis'^^s'} = {akis,akis'} =0 
|akis,a^s/| = <5kiki<5ss'- 



(11) 
(12) 



The distinction between BCS pairs and CPs holds for the original or "ordinary" CPs [H] as in ^ and 

It also applies to the generalized BS CPs defined more consistently without excluding two-hole pairs when 

the lower limit in ([1]) is taken as \/kp — kf-^ as in BCS theory |T7]. 



The BCS-pair annihilation and creation operators for any K > can be defined quite generally as 



OkK = OkaXakiT 



and 



\K — "ki 



(13) 



where aj^.^^ and Ok^s obey (fTT|) and (fT2|) . and as before k = ^(ki — and K = ki +k2 are the relative and 
total (or center-of-mass) momentum wavevectors, respectively, associated with two fermions with wavevectors 



ki K/2 + k 



and 



ka ==K/2-k. 



(14) 



Using the same techniques to derive © to ([5]) valid for K = 0, the operators 6kK and 6j^j^ are found to 
satisfy: a) the "pseudo-boson" commutation relations 



6kK, 6jj.,K — (1 - riK/2-ki ~ "K/2-|-kT)'5kk' 



(15) 



"kK' "k'K 



[6kK,6k'K] = 



(16) 



where n. 



K/2±k3 

commutation relation 



'^K/2±ks'^K/2±ks are fermion number operators; as well as b) a "pseudo-fermion" anti- 



{6kK,6k'K} = 26kK6k'K(l - ^kk')- 



(17) 



Our only restriction was that K = ki + ka = k'j^ + kj. If K = so that ki — — ka — k (the only case 
considered by BCS), and caUing 6kK=o = 6k, etc., (fT5|l to p7|) become ([6]) to ([8]), as they should. So, neither 
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BCS pairs with K > are bosons as the relation (fTS]) contains additional terms not [51] present in the the 
usual boson commutation relations analogous to ([9]). 

To our knowledge, no one has yet succeeded in constructing CP creation and annihilation operators that 
obey Bose commutation relations, starting from fermion creation aj^^^ and annihilation ay^^g operators, as 
is postulated in Refs. [U [3| in a generalized EEC theory that contains BCS theory as a special case. This 
postulation is grounded in magnetic-flux quantization experiments [S^-IMj that establish the presence of 
charged pairs — albeit without being able [55] to specify the sign of those charges. However, although the 
eigenvalues of ^^x^kK are or 1 in keeping with the Pauli Exclusion Principle, those of X^k^kK^kK are 
evidently 0, 1, 2, ... because of the indefinitely many values taken on by the summation index k. This implies 
BE statistics and corroborates the qualitative conclusions reached above. A discussion in greater detail of 
this is found in Refs. [56] [57j. 



4 Conclusion 



In conclusion, in the thermodynamic limit any number of CPs with a definite total (or center-of-mass) 
wavevector K can occupy a single state of CP energy £k and thus obey BE statistics. Hence, CPs can 
undergo a BEC. All this holds for any coupling, i.e., regardless of the size of the CPs and of their mutual 
overlap. As in liquid ^He, the BEC statistical mechanism may serve as a starting point — even before the 
effects of interboson interactions (excluded also in BCS theory) are accounted [58j for — in constructing a 
microscopic theory of superconductivity, or of neutral-fcrmion superfluidity, that will eventually lead to Tc 
values calculated from first principles without adjustable parameters and that hopefully agree with observed 
ones. 
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